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Abstract. The importance of a node in a directed graph can be measured by its PageRank. 
The PageRank of a node is used in a number of application contexts - including web-page ranking — 
and can be interpreted as the average portion of time spent at the node by an infinite random walk 
in the graph. We consider the problem of maximizing the PageRank of a node by selecting some of 
the edges from a set of edges that are under our control. By applying results from Markov decision 
theory, we show that an optimal solution to this problem can be found in polynomial time. Our 
approach provides a linear programming formulation that can then be solved by standard techniques. 
We also show in the paper that, under the slight modification of the problem for which we are given 
mutually exclusive pairs of edges, the problem of PageRank optimization becomes NP-hard. 
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1. Introduction. The importance of a node in a directed graph can be mea- 
sured by its PageRank. The PageRank of a node can be interpreted as the average 
portion of time spent at the node by an infinite random walk [14j . PageRank is tradi- 
tionally applied for ordering web-pages of the web-graph in which the edges represent 
hyperlinks, but it also has many other applications [2], for example, in bibliometrics, 
ecosystems, spam detection, web-crawling, semantic networks, keyword searching in 
relational databases and word dependencies in natural language processing. 

It is of natural interest to search for the maximum or minimum PageRank that 
a node (web-page) can have depending on the presence or absence of some of the 
edges (hyperlinks) in the graph. For example, since PageRank is used for ordering 
web search results, a web-master could be interested in increasing the PageRank of 
some of his web-pages by suitably placing hyperlinks on his own site or by buying 
advertisements or making alliances with other sites [TJ [5] • Another motivation is that 
of estimating the PageRank of a node in the presence of missing information on the 
graph structure. If some of the links on the internet are broken, for example, because 
the server is down or there are network traffic problems, we may have only partial 
information on the link structure of the web-graph. However, we may still want to 
estimate the PageRank of a web-page by computing the maximum and minimum 
PageRank that the node may possibly have depending on the presence or absence of 
the hidden hyperlinks [12] . These hidden edges are often referred to as "fragile links" . 

It is known that if we place a new edge in a directed graph, then the PageRank of 
the target node can only increase. Optimal linkage strategies are known for the case 
in which we want to optimize the PageRank of a node and we only have access to the 
edges starting from this node pQ. This first result has later been generalized to the 
case for which we are allowed to configure all of the edges starting from a given set 
of nodes [8] . The general problem of optimizing the PageRank of a node in the case 
where we are allowed to decide the absence or presence of the edges in a given subset 
of edges is proposed by Ishii and Tempo |12j . They are motivated by the problem 
of "fragile links" and they mention the lack of polynomial time algorithms to this 
problem. Then, using interval matrices, they propose an approximate solution |12j . 
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In the paper we show that this optimization problem can be efficiently formu- 
lated as a Markov decision process (MDP), more precisely, as a stochastic shortest 
path (SSP) problem, and that it can therefore be solved in polynomial time. Our 
proof of polynomial time computability is based on a linear programming formula- 
tion, but we also propose an iterative algorithm that has appealing properties, for 
example, it always finds an optimal solution in finite steps and, under some assump- 
tions, its running time is polynomial. In a last section, we prove that under the slight 
modification of the problem for which we are given mutually exclusive constraints 
between pairs of edges, the problem of PageRank optimization becomes NP-hard. 

2. Definitions and Preliminaries. In this section we first give on overview on 
PageRank computation, then we define the problem of optimizing the PageRank of a 
given node via a subset of edges, called fragile links. We also give a brief introduction 
to stochastic shortest path problems, since they constitute the basis of our solution. 

2.1. PageRank. Let Q = (V, £) be a directed graph, where V = {1, . . . , n} is 

the set of vertices and £ C V x V is the set of edges. First, for simplicity, we assume 
that Q is strongly connected. The adjacency matrix of Q is denoted by A. Since Q is 
strongly connected, A is irreducible. We are going to consider a random walk on the 
graph defined as follows. If we are in node i, in the next step we will go to node j 
with probability l/deg(i) if j is an out-neighbor of i, where deg(-) denotes out-degrcc. 
This defines a Markov chain on the nodes of the graph with transition-matrix 

P=(D A 1 A) T with Da = diag(At) (2.1) 

where 1 = (1, . . . , 1} T is the all-one vector and diag(-) is an operator that creates a 
diagonal matrix from a vector, more precisely, (Da)u = {At)i = deg(i). Note that 
P is a column (left) stochastic matrix and the chain can be interpreted as an infinite 
uniform random walk on the graph, or in the network terminology, a random surfing. 

The PageRank vector, 7r, of the graph is defined as the stationary (steady-state) 
distribution of the above described homogeneous Markov chain, more precisely, as 

Ptv = tv, (2.2) 

where 7r is non-negative and 7r T l = 1. Since P is an irreducible stochastic matrix, 
we know, e.g., from the Perron- Frobenius theorem, that 7r exists and it is unique. 

Now, we turn to the more general case, when we do not assume that Q is strongly 
connected, it can be an arbitrary directed graph. In this case, there may be nodes 
which do not have any outgoing edges. They are usually referred to as "dangling" 
nodes. There are many ways to handle dangling nodes [2J, e.g., we can delete them, 
we can add a self-loop to them, each dangling node can be linked to an ideal node 
(sink) or we can connect each dangling node to every other node. This last solution 
can be interpreted as restarting the walk from a random starting state once we have 
reached a dangling node. Henceforth, we will assume that we have already dealt with 
the dangling nodes and, therefore, every node has at least one outgoing edge. 

We can define a Markov chain similarly to (|2.1|) . but this chain may not have a 
unique stationary distribution. The solve this problem, the PageRank vector, tt, of Q 
is defined as the stationary distribution of the "Google matrix" [T3] defined as 

G = (l-c)P + czt T , (2.3) 

where z is a positive "personalization" vector satisfying z T l = 1, and c 6 (0, 1) is a 
"damping" constant. In practice, values between 0.1 and 0.15 are usually applied for 
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c and z = (l/n)l [2]. The Markov chain defined by G is regular (or crgodic) that 
is irreducible and aperiodic, therefore, its stationary distribution, tt, uniquely exists 
and, moreover, the Markov chain converges to tt from any initial distribution [15j . 

The idea of PageRank is that Tz(i) can be interpreted as the "importance" of i. 
Thus, tt defines a linear order on the nodes by treating i < j if and only if tt(i) < tt(J). 

The PageRank vector can be approximated by the power method, namely, by 

X n +1 = Gx n , (2.4) 

where xq is an arbitrary stochastic vector. It can also be directly computed by 

tt = c {I -{l-c)P)- 1 z 1 (2.5) 

where I denotes an n x n identity matrix. Since c G (0,1) and P is stochastic, it follows 
that matrix I — (1 — c)P is strictly diagonally dominant, therefore, it is invcrtiblc. 

2.2. PageRank Optimization. We will investigate the problem when a subset 
of links are "fragile" , we do not know their exact values or we have control over them, 
and we want to compute the maximum (or minimum) PageRank that a specific node 
can have [12]. More precisely, we are given a digraph Q = (V, £), a node v £ V and a 
set T C £ corresponding to those edges which are under our control. It means that we 
can choose which edges in T are present and which are absent, but the edges in £ \ T 
are fixed. We will call any T+ C T a configuration of fragile links: T + determines 
those edges that we add to the graph, while T- = T \ J-+ denotes those edges which 
we remove. The PageRank of node v under the J-+ configuration is defined as the 
PageRank of v with respect to the graph Qo = (V,£\J--). The problem is that 
how should we configure the fragile links, in order to maximize (or minimize) the 
PageRank of vl The Max-PagcRank problem can be summarized as follows: 



The Max-PageRank Problem 

Instance: A digraph Q = (V,£), a node v £ V and a set of controllable edges T C £. 
Optional: A damping constant c G (0, 1) and a stochastic personalization vector z. 
Task: Compute the maximum possible PageRank of v by changing the edges in T 

and provide a configuration of edges in T for which the maximum is taken. 



The Min-PagcRank problem can be stated similarly. We will concentrate on the 
Max problem but a straightforward modification of our approach can deal with the 
Min problem, as well. We will show that Max-PageRank can be solved in polynomial 
time even if the damping constant and the personalization vector are part of the input . 

There are finitely many configurations, consequently, we can try to compute them 
one-by-one. If we have d fragile links, then there are 2 d possible graphs that we can 
achieve. Given a graph, we can compute its PageRank vector in 0(n 3 ), via a matrix 
inversion^. The resulting "brute force" algorithm has 0(n 3 2 d ) time complexity. 

We note that if the graph was undirected, the Max-PageRank problem would be 
easy. We know [17] that a random walk on an undirected graph (a time-reversible 
Markov chain) has the stationary distribution Tr(i) = deg(i)/2m for all i, where m 
denotes the number of edges and deg(i) is the degree of i. Hence, in order to maximize 
the PageRank of v, we should keep edge £ J 7 if and only if i — v or j = v. 



1 It can be done a little bit faster, in 0(n 2 376 ), using the Coppersmith- Winograd algorithm [7]. 
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2.3. Stochastic Shortest Path Problems. In this section we give an overview 
on stochastic shortest path problems, since our solutions arc built upon their theory. 

Stochastic shortest path (SSP) problems are generalizations of (deterministic) 
shortest path problems [U[S]. In an SSP problem the transitions between the nodes 
are uncertain, but we have some control over their probability distribution. We aim 
at finding a control policy (a function from nodes to controls) such that by applying 
this policy we can reach a given target state with probability one while minimizing 
the expected sum of costs (which depends on the applied controls), as well. SSP prob- 
lems are special undiscounted Markov decision processes (MDPs) with finite state and 
action spaces and with an absorbing, cost-free terminal state. They are of high prac- 
tical importance since they arise in many real-world domains, for example, several 
problems in operations research can be formulated as an SSP problem or contains it 
as a sub-problem. Typical examples are routing and resource allocation problems. 

2.3.1. Problem Formulation. An SSP problem can be stated as follows. We 
have given a finite set of states, S, and a finite set of control actions, U. For simplicity 
we assume that § = {1, . . . , n, n + 1}, where r = n + 1 is a special state, the target or 
termination state. In each state i we can choose an action u G U(i), where U(i) CD 
is the set of allowed actions in state i. After the action was chosen, the system moves 
to state j with probability p(j \ i,u) and we incur cost g(i,u,j). The cost function is 
real valued and the transition-probabilities are, of course, nonnegative as well as they 
sum to one. Namely, for all state i, j and action u £ U(i), p(j \i,u) > and 

n+l 

£p(7|i,«) = l. (2.6) 
j'=i 

The target state is absorbing and cost-free that is, if we reach state t, we remain there 
forever without incurring any more costs. More precisely, for all action u G U(j), 

P (t\t,u) = 1, 9 {t,u,t) = 0. (2.7) 

The problem is to find a control policy such that it reaches state r with probability 
one and minimizes the expected cumulative costs, as well. A (stationary, Markovian) 
deterministic policy is a function from states to actions: /i : § — > U. A randomized 
policy can be formulated as fi : § — > A(U), where A(U) denotes the set of all probabil- 
ity distributions over set U. It can be shown that every policy induces a Markov chain 
on the state space [9] . A policy is called proper if, using this policy, the termination 
state will be reached with probability one, and it is improper otherwise. The value or 
cost-to-go function of policy /i gives us the expected total costs that we incur starting 
from a state and following fi thereafter; that is, for all state i, it is defined as 



J^(i) = Urn 

k — >oo 



fc-1 

^2g(i t ,u t ,it+i) 



(2.8) 



where i t and u t are random variables representing the state and the action taken at 
time t, respectively. Naturally, state it+i is of distribution p(- \ it,Ut) and action u t is 
of distribution fJ-(it), or u t = n{it) if the applied policy was deterministic. 

Applying a proper policy, we arrive at a finite horizon problem, however, the 
length of the horizon may be random and may depend on the applied policy, as well. 
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2.3.2. Bellman Equations. Wc say that policy /Xi is better than or equal to 
policy H2 if for all state i, J Ml (i) < J^ 2 (i). A policy is (uniformly) optimal if it is 
better than or equal to all other policies. There may be many optimal policies, but 
assuming that (Al) there exists at least one proper policy and (A2) every improper 
policy yields infinite cost for at least one initial state, they all share the same unique 
optimal value function, J*. Then, function J* is the unique solution of the Bellman 
equation, TJ* = J*, where T is the Bellman operator [5]; that is, for all i, 

n+1 

(TJ)(i)= imn p{j\i,u)\g{i,u,j) + J(j)]. (2.9) 

u&AU,) — ■ L J 

3=1 

The Bellman operator of a (randomized) policy /i is defined for all state i as 

n+1 

(T a J)(i)= ^2n(i,u)^2pU\i,u)[g(i,u,j) + J(j)], (2.10) 
ueu(i) j=i 

where fj,(i,u) is the probability that policy \i chooses action u in state i. Naturally, 
this definition is valid for deterministic policies, as well, since they are special cases 
of randomized ones. Given the assumptions above we also have T u J M = J M and the 
standard value iteration converges in SSPs [1] for all initial vector J, that is 

lim T*J = J fl , limT k J = J*. (2.11) 

k — >oo k — >oc 

It is also known that a policy [i is optimal if and only if T U J* = TJ* (= J*). 
Moreover, the Bellman operator, T, and the Bellman operator of any policy fi, T u , 
are monotone and, assuming that (APP) all allowed control policies are proper, T 
and T u are contractions with respect to a weighted maximum norm [5]. 

From a given value function J, it is straightforward to get a policy, for example, 
by applying a greedy (with respect to J) policy; that is, for all state i, 

n+1 

H(i) 6 argmin V p(j | i,u) g(i,u,j) + J(j) . (2.12) 
uew(i) j=1 J 

The policy iteration (PI) algorithm works as follows. We start with an arbitrary 
proper policy, fiQ . In iteration k we first evaluate the actual policy, /Lt^ , by solving the 
linear system, T Hk J tlk = J flk , and then we improve the policy by defining fik+i as the 
greedy policy with respect to J Mfe . The algorithm terminates if J^ h = J AIfc+1 . It is 
known that, assuming (Al) and (A2), PI generates an improving sequence of proper 
policies and finds an optimal solution in finite number of iterations [5] . 

There are many solution methods for solving MDPs, e.g., in the fields of reinforce- 
ment learning and dynamic programming. Most of these algorithms aim at finding 
(or approximating) the optimal value function, since it immediately gives an optimal 
policy, as well. Such methods are, e.g., value iteration, policy iteration, Gauss-Seidel 
method, Q-learning, Q(A), SARSA and TD(A): temporal difference learning [5l l9l 120] . 

2.3.3. Computational Complexity. It is known that all of the three classical 
variants of MDPs (finite horizon, infinite horizon discounted cost and infinite horizon 
average cost) can be solved in polynomial time [TB]. Moreover, these classes of prob- 
lems are P-complete fT5]. In the case of SSP problems, it is know that they can be 
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reformulated as linear programming (LP) problems [5]. More precisely, the optimal 
cost-to-go function, J*(l), . . . , J*(n), solves the following LP in variables X\, . . . , x n : 



maximize 



Xi (2.13a) 

i=l 

n+1 

subject to Xi < p(j I i, u) g(i,u,j) + Xj (2.13b) 

for all state i and for all action u S Note that x n +i is fixed at zero. This LP 

has n variables and 0(nm) constraints, where m is the maximum number of allowed 
actions per state. Knowing that an LP can be solved in polynomial time (in the 
number of variables, the number of constraints and the binary size of the input), this 
reformulation already provides a way to solve an SSP problem in polynomial time. 

Though, value iteration usually requires an infinite number of steps to obtain the 
optimal value function, stopping rules can be defined which guarantee that the greedy 
policy with respect to the current value function approximation will be optimal [16| . 

Assuming that all policies are proper (APP), the state space can be partitioned 
into nonempty subsets S\, . . . , S r such that for any 1 < q < r, state i £ S q and action 
u G U(i), there exists some j <E {t} U Si U • • ■ U S q -i such that p{j \ i,u) > 0. Then, 
if assumption (APP) holds, value iteration can find an optimal policy after a number 
of iterations that is bounded by a polynomial in L (the binary input size) and r\~ 2r , 
where r\ is the smallest positive transition probability of the MDP [21] . 

Since policy iteration converges no more slowly than value iteration |19| . it also 
terminates in iterations bounded by a polynomial in L and rj~ 2r , assuming (APP). 



3. Optimizing PageRank as a Markov Decision Process. Before we prove 
that efficient (polynomial time) algorithms to the Max-PageRank problem do exist, 
first, we recall a basic fact about stationary distributions of Markov chains. 

Let (Xq, Xi,...) denote a homogeneous Markov chain defined on a finite set f£. 
The expected first return time of a state i € f2 is defined as follows 

<p(i) = E [ min { t > 1 : X t = i } | X Q = i } . (3.1) 

If state i is recurrent, then tp(i) is finite. Moreover, if the chain is irreducible, 

for all state i, where 7r is the stationary distribution of the chain |15j . This naturally 
generalizes to unichain processes, viz., when we have a single communicating class 
of states and possibly some transient states. In this case we need the convention 
that l/oo = 0, since the expected first return time to transient states is oo. Hence, 
the stationary distribution of state i can be interpreted as the average portion of 
time spent in i during an infinite random walk. It follows from equation (|3.2[) that 
maximizing (minimizing) the PageRank of a given node is equivalent to the problem 
of minimizing (maximizing) the expected first return time to this node. 

We are going to show that the Max-PageRank problem can be efficiently formu- 
lated as a Markov decision process (MDP), more precisely, as a stochastic shortest 
path (SSP) problem [H[5], where "efficiently" means that in polynomial time. First, 
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we will consider the PageRank optimization without damping, namely c = 0, but later, 
we will extend the model to the case of damping and personalization, as well. 

We will start with a simple, but intuitive reformulation of the problem. Though 
this reformulation will not ensure that Max-PageRank can be solved in polynomial 
time, it is good to demonstrate the main ideas and to motivate the refined solution. 

3.1. Assumptions. First, we will make two assumptions, in order to simplify 
to presentation, but later, in the main theorem, they will be relaxed. 

(AD) Dangling Nodes Assumption: We assume that there is a fixed (not fragile) 
outgoing edge from each node of the graph. This assumption guarantees that 
there are no dangling nodes and there are no nodes with only fragile links. 

(AR) Reachability Assumption: We also assume that for at least one configuration 
of fragile links we have a unichain process and node v is recurrent, namely, we 
can reach node v with positive probability from all nodes. This assumption 
is required to have a well-defined PageRank for at least one configuration. In 
our SSP formulation this assumption will be equivalent to assuming that there 
is at least one proper policy, which statement can be checked in polynomial 
time [S]. In case of damping this assumption is automatically satisfied, since 
then the Markov chain is irreducible, and hence unichain, no matter how we 
configure the fragile links, which implies that all policies are proper. 

3.2. Simple SSP Formulation. First, let us consider an instance of the Max- 
PageRank Problem. We will build an associated SSP problem as follows. The states 
of the MDP are the nodes of the graph, except for v which we "split" into two parts 
and replace by two new states: v s and Vt- Intuitively, state v s will be our "starting" 
state: it has all the outgoing edges of v (both fixed and fragile), but it does not have 
any incoming edges. The "target" state will be vt: it has all the incoming edges of 
node v and, additionally, it has only one outgoing edge: a self- loop. Note that r = Vt, 
namely, vt is the absorbing termination state of the SSP problem. 

An action in state i is to select a subset of fragile links (starting from i) which 
we "turn on" (activate). All other fragile links from this state will be "turned off" 
(deactivated). Therefore, for each state i, the allowed set of actions is U{i) = V{Ti), 
where V denotes the power set and Ti is the set of all outgoing fragile edges from i. 

Let us assume that we are in state i, where there are > 1 fixed outgoing edges 
and we have activated bi(u) > fragile links, determined by action u € U(i). Then, 
the transition- probability to all state j that can be reached from state i using a fixed 
or an activated fragile link (an out-neighbor of i w.r.t. the current configuration) is 



We define the immediate-cost of all control actions as one, except for the actions 
taken at the cost-free target state. Therefore, the immediate-cost function is 



for all state i, j and action u. Note that taking an action can be interpreted as 
performing a step in the random walk. Therefore, the expected cumulative cost of 
starting from state v s until we reach the target state Vt is equal to the expected 




(3.4) 
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number of steps until we first return to node v according to our original random walk. 
It follows, that the above defined SSP formalizes the problem of minimizing (via an 
appropriate configuration of fragile links) the expected first return time to state v. 
Consequently, its solution is equivalent to maximizing the PageRank of node v. 




Each allowed deterministic policy [i defines a potential way to configure the fragile 
links. Moreover, the v s component of the cost-to-go function, J M (w s ), is the expected 
first return time to v using the fragile link configuration of fi. Therefore, we can 
compute the PageRank of node v, assuming the fragile link configuration of /i, as 

J^(Vs) 

where we applied the convention of l/oo = 0, which is needed when v is not recurrent 
under the configuration of \i. Thus, the maximal PageRank of v is 1/J*(v s ). 

Most solution algorithms compute the optimal cost-to-go function, J*, but even 
if we use a direct policy search method, it is still easy to get back the value function 
of the policy. We can compute, e.g., the expected first return time if we configure the 
fragile links according to policy \i as follows. For simplicity, assume that v s = 1, then 

J"(l) = l T (J-P M )- 1 e 1 , (3.6) 

where ej is j-th canonical basis vector, I is an n x n identity matrix and P M is the 
substochastic transition matrix of the SSP problem without the row and column of the 
target state, Vt, if we configure the fragile links according to policy \i. Note that 

oo 
n=0 

and we know that this Neumann series converges if g(P fl ) < 1, where g(-) denotes the 
spectral radius. Therefore, (/ — P^) _1 is well-defined for all proper policies, since it is 
easy to see that the properness of policy fi is equivalent to g{P^) < 1. 

3.3. Computational Complexity. It is known that MDPs can be solved in 
polynomial time in the number of states, N, and the maximum number of actions per 
state, M (and the maximum number of bits required to represent the components, 
L), e.g., by linear programming [TBI HB]- The size of the state space of the current 
formulation is N = n + 1, where n is the number of vertices of the original graph, 
but, unfortunately, its action space does not have a polynomial size. For example, 
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if we have maximum m fragile links leaving a node, we had 2 m possible actions to 
take, namely, we could switch each fragile link independently on or off, consequently, 
M = 2 m . Since m = O(n), from the current reformulation of problem, we have that 
there is a solution which is polynomial but in 2™, which is obviously not good enough. 
However, we can notice that if we restrict the maximum number of fragile links per 
node to a constant, k, then we could have a solution which is polynomial in n (since 
the maximum number of actions per state becomes constant: 2 k ). This motivates our 
refined solution, in which we are going to reduce the maximum number of allowed 
actions per state to two while only slightly increasing the number of states. 

3.4. Refined SSP Formulation. We are going to modify our previous SSP 
formulation. The key idea will be to introduce an auxiliary state for each fragile link. 
Therefore, if we had a fragile link from node i to node j in the original graph, we 
place an artificial state, fy, "between" them in the refined reformulation. The refined 
transition-probabilities are as follows. Let us assume that in node i there were a, > 1 
fixed outgoing edges and bi > fragile links. Now, in the refined formulation, in 
state i we have only one available action which brings us uniformly, with l/(a; + 6,) 
probability, to state j or to state /y depending respectively on whether there was a 
fixed or a fragile link between i and j. Notice that this probability is independent of 
how many fragile links are turned on, it is always the same. In each auxiliary state 
fij we have two possible actions: we could either turn the fragile link on or off. If our 
action was "on" then we went with probability one to state j, however, if our action 
was "off", then we went back with probability one to state i (Figure l3.2p . 




Figure 3.2. An example for inserting auxiliary sates for fragile links. The left hand side 
presents the original situation, in which dashed edges are fragile links. The right hand side shows the 
refined reformulation, where the dotted edges represent possible actions with deterministic outcomes. 

We should check whether the transition-probabilities between the original nodes 
of graph are not affected by this reformulation. Suppose, we are in node i, where 
there are a fixed and b fragile link^, and we have turned k of the fragile links on. 
Then, the transition-probability to each node j, which can be reached via a fixed or an 
activated fragile link, should be l/(a + k). In our refined reformulation, the immediate 
transition-probability from state i to state j is l/(a+b), however, we should not forget 
about those b — k auxiliary nodes in which the fragile links are deactivated and which 
lead back to state i with probability one, since, after we have returned to state i we 



2 For simplicity, now we do not denote their dependence on node i. 
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have again l/(a + b) probability to go to state j an so on. Now, wc arc going compute 
the probability of eventually arriving at state j if we start in state i. 

To simplify the calculations, let us temporarily replace each edge leading to an 
auxiliary sate corresponding to a deactivated fragile link with a self-loop. We can 
safely do so, since these states lead back to state i with probability one, therefore, the 
probability of eventually arriving at j does not change by this. Then, we have 

¥(3t:X t =j\Vs<t:X s = i) = (3.8a) 
= ^F{Xt=j\X t -i = i) l[¥(X s = i\X s _ 1 = i) = (3.8b) 

t=l s=l 




(3.8c) 



With this, we have proved that the probability of eventually arriving at state j if 
we start in state i, before arriving at any (non-auxiliary) node / that was reachable 
via a fixed or a fragile link from i in the original graph, is the same as the one-step 
transition-probability was from i to j according to the original random walk. 

However, we should be careful, since we might have performed several steps in 
the auxiliary nodes before we finally arrived at state j. Fortunately, this phenomenon 
does not ruin our ability to optimize the expected first return time to state v in the 
original graph, since we count the steps with the help of the immediate-cost function, 
which can be refined according to our needs. All we have to do is to allocate zero cost 
to those actions which lead us to auxiliary states. More precisely, 

/. -n • / if i = v t or j = f a or u = "off", 
S(w) = (l otherwise, (3 ' 9) 

for all state i, j, I and action u. Consequently, we only incur cost if we directly go 
from state i to state j, without visiting an auxiliary node (viz., it was a fixed link), 
or if we go to state j via an activated fragile link, since we have g(fij,u,j) = 1 if 
u = "on" . It is easy to see that in this way we only count the steps of the original 
random walk and, e.g., it does not matter how many times do we visit auxiliary nodes. 

This reformulation also has the nice property that J^fos) is the expected first 
return time to node v in the original random walk, in case we have configured the 
fragile links according to policy fj.. The minimum expected first return time that we 
can achieve with suitably setting the fragile links is J*(v s ), where J* is the optimal 
cost-to-go function, consequently, the maximum PageRank of node v is 1/ J*(v s ). 

It is also easy to see that if we want to compute the minimum possible PageRank 
of node v, then we should simply define a new immediate-cost function as g = —g, 
where g is defined by equation (|3.9p . If the optimal cost-to-go function of this modified 
SSP problem is J*, then the minimum PageRank that node v can have is l/\ J*(v s )\. 

3.5. Computational Complexity. The number of states of this formulation 
isN = n + d-\-l, where n is the number of nodes of the original graph and d is the 
number of fragile links. Moreover, the maximum number of allowed actions per state 
is M = 2, therefore, this SSP formulation provides a proof that, assuming (AD) and 
(AR), Max-PagcRank can be solved in polynomial time in the size of the problem. 
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The resulted SSP problem can be reformulated as a linear program, namely, the 
optimal cost-to-go function solves the following LP in variables Xi and , 



maximize 



subject to 



E x * + E 

Xj-j ^ Xi , and 



Xjj < Xj + 1 , and 



Xi < 



1 

deg(i) 



E & + + E 



(3.10a) 
(3.10b) 
(3.10c) 



for alH G V \ and (i, j) G J 7 , where Xi is the cost-to-go of state i, Xij relates to 
the auxiliary states of the fragile edges, and deg(-) denotes out-degree including both 
fixed and fragile links (independently of the configuration). Note that we can only 
apply this LP after state v was "splitted" into a starting and a target state, hence, 
v Sl v t G V. Also note that the value of the target state, x Vt , is fixed at zero. 

3.6. Handling Dangling Nodes. Now, we will show that assumption (AD) 
can be omitted and our result is independent of how the dangling nodes are handled. 

Suppose that we have chosen a rule according to which the dangling nodes are 
handled, e.g., as discussed by Berkhin [2]. Then, in case (AD) is not satisfied, we 
can simply apply this rule to the dangling nodes before the optimization. However, 
we may still have problems with the nodes which only have fragile links, since they 
are latent dangling nodes, namely, they become dangling nodes if we deactivate all of 
their outgoing edges. We call them "fragile nodes" . Notice that we can safely restrict 
the optimization in a way that maximum one of the links can be activated in each 
fragile node. This does not affect the optimal PageRank, since the only one allowed 
link should point to a node that has the smallest expected hitting time to v. Even if 
there are several nodes with the same value, we can select one of them arbitrarily. 

It may also be the case that deactivating all of the edges is the optimal solution, 
e.g., if the fragile links lead to nodes that have very large hitting times to v. In this 
case, we should have an action that has the same effect as the dangling node handling 
rule. Consequently, in case we have a fragile node that has m fragile links, we will 
have in + 1 available actions: u±, . . . , m to +i- If Uj is selected, where 1 < j < m, it 
means that only the j-th fragile link is activated and all other links are deactivated, 
while if u m+ i is selected, it means that all of the fragile links are deactivated and 
auxiliary links are introduced according to the selected dangling node handling rule. 

If we treat the fragile nodes this way, we still arrive at an MDP which has states 
and actions polynomial in n and d, therefore, Max-PagcRank can be solved in poly- 
nomial time even if (AD) is not satisfied and independently of the applied rule. The 
modification of the LP formulation if fragile nodes are allowed is straightforward. 

3.7. Damping and Personalization. Now, we are going to extend our refined 
SSP formulation, in order to handle damping, as well. For the sake of simplicity, we 
will assume (AD), but it is easy remove it in a similar way as it was presented in 
Section 13.61 Note that assumption ( AR) is always satisfied in case of damping. 

Damping can be interpreted as in each step we continue the random walk with 
probability 1 — c and we restart it with probability c, where c G (0, 1) is a damp- 
ing constant. In this latter case, we choose the new starting state according to the 
probability distribution of a positive and stochastic personalization vector z. In order 
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to model this, we introduce a new global auxiliary state, q, which we will call the 
"teleportation" state, since random restarting is sometimes called telcportation |14j . 



fragile link 




Figure 3.3. An illustration of damping: the substructure of a node of the original digraph. 
Circles represent states and boxes represent actions. The double circles indicate that in fact the 
actions do not directly lead to these states, but to their similarly defined substructures. State q 
denotes the global "teleportation" state from which the states can be reached with distribution defined 
by the personalization vector. Dashed edges help determining zero cost events: if a state-action-state 
path has only dashed edges, it means that this triple has zero cost, otherwise, its cost is one. 



In order to take the effect of damping into account in each step, we place a new 
auxiliary state hi "before" each (non-auxiliary) state i (see Figure [3T3]) . Each action 
that leaded to i in the previous formulation now leads to hi. In hi we have only one 
available action ( "nop" abbreviating "no operation" ) which brings us to node i with 
probability 1 — c and to the teleportation state, q, with probability c, except for the 
target state, Vt, for which h Vt leads with probability one to v t . In q, we have one 
available action which brings us with z distribution to the newly defined nodes, 

!z(i) if i ^ v s and i ^ v t 
z(v) if i = v t (3-11) 
if i = v s . 

All other transition-probabilities from q are zero. Regarding the immediate-cost 
function: it is easy to see that we should not count the steps when we move through 
hi, therefore, g(hi,u,i) = and g(hi,u,q) = 0. However, we should count when we 
move out from the teleportation state, i.e., g(q,u,i) = 1 for all state i and action u. 

3.8. Computational Complexity. In the current variant, which also takes 
damping into account, the size of the state space is N = 2n + d + 2, and we still have 
maximum 2 actions per state, therefore, it can also be solved in polynomial time. 
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In this case, the LP formulation of finding the optimal cost-to-go is 



maximize 



x t + Xi) + 2J 3 



(3.12a) 



subject to 




and Xi < (1 — c) Xi + cx, 



(3.12b) 



and x q < Zi (£i + 1) , 



(3.12c) 



1 



E + + E 



(3.12d) 



Xi < 



deg(i) 



for all i G V \ {vt} and G J 7 , where Zi = p(hi \ q), Xi denotes the cost-to-go of 
state hi and x q is the value of the teleportation state, q. All other notation are the 
same as in LP (|3.10[) and we also have that the values x Vt and fixed at zero. 

We arrived at an LP problem with 0{n + d) variables and 0(n + d) constraints. 
Given an LP with k variables and 0{k) constraints, it can be solved in 0(k 3 L), where 
L is the binary input size (for rational coefficients) or in 0(fc 3 logi), where £ is the 
desired precision [TT]. Therefore, Max-PageRank can be solved using 0((n + d) 3 L) 
operations under the Turing model of computation. We can conclude that 

Theorem 3.1 (Main Theorem). The Max-PageRank Problem can be solved 
in polynomial time even if the damping constant and the personalization vector are 
part of the input and independently of the way the dangling nodes are handled. 

Note that assumptions (AD) and (AR) are not needed for this theorem, since dan- 
gling and fragile nodes can be treated as discussed in Section 13.61 (without increasing 
the complexity of the problem) and, in case of damping, all policies are proper. 

Assuming the damping constant, c, and the personalization vector, z, can be 
represented using a number of bits polynomial in n, which is the case in practice, 
since c is usually 0.1 or 0.15 and z = (1/n) 1 [2], we arrive at a strongly polynomial 
time solution, because the other coefficients can be represented using O(logn) bits. 

3.9. State Space Reduction. In the last SSP formulation in 2n + 1 states 
there is no real choice in (there is only one available action), therefore, the state space 
may be reduced. In this complementary section we are going to show that given an 
SSP problem with N = r + s states, in which in r states there is only one available 
action, we can "safely" reduce the number of states to s. More precisely, we will prove 
that we can construct another SSP problem with only s states which is "compatible" 
with the original one in the sense that there is a one-to-one correspondence between 
the policies of the reduced and the original problems, and the value functions of the 
policies (restricted to the remaining s states) are the same in both problems. Hence, 
finding an optimal policy for the reduced problem is equivalent to solving the original 
one. The computational complexity of the construction is 0(r 3 +r 2 sm+ s 2 rm), where 
m denotes the maximum number of allowed actions per state. 

We will apply immediate-cost functions of the form g : § x U — > R. If we had a 
cost function that also depended on the arrival state, we could redefine it by 



N 
3 = 1 



(3.13) 
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which would not affect the outcome of the optimization [3] . Note that computing the 



new cost function can be obtained using 0(N m) = 0(r m- 



- s m) operations. 



We will call the states in which there is only one available control action as 
"non-decision" states, while the other states will be called "decision" states. By 
convention, we classify the target state, r, as a decision state. We will assume, w.l.o.g., 
that the (indices of the) non-decision states are 1, . . . , r and the decision states are 
?* + 1, ...,?* + s. Finally, we will assume that there exists at least one proper policy. 

3.9.1. Construction. Notice that the transition-matrix of the Markov chain 
induced by (any) control policy [i of the original SSP problem looks like 



Qo Qii 



(3.14) 



where Ro £ W xr describes the transition-probabilities between the non-decision 
states; Qo £ M. sxr contains the transitions from the non-decision states to the de- 
cision states; Q^ £ R sxs describes the transitions between the decision states and, 
finally, i? M £ M. rxs contains the transitions from the decision states to the non-decision 
states. Notice that Ro and Qo do not depend on the applied control policy. 

In the reduced problem we will only keep the s decision states and remove the r 
non-decision states. We are going to redefine the transition-probabilities (and later 
the costs) between the decision states as if we would "simulate" the progress of the 
system through the non-decision states until we finally arrive at a decision state. In 
order to calculate the probabilities of arriving at specific decision states if we started 
in specific non-decision states, we can define a new Markov chain as 



Po = 



Ro 
Qo I 



(3.15) 



where is a r x s zero matrix and I is an s x s identity matrix. We can interpret this 
matrix as if we would replace each decision state by an absorbing state. There is at 
least one proper policy, Ro and Qo are the same for all policies and the target state 
is a decision state, therefore, Rq converges to the zero matrix as k — > oo, thus 



lim P k = 




Q* I 



(3.16) 



where Q* contains the arrival distributions to the decision states if we started in one of 
the non-decision states. More precisely, Q*- is the probability of arriving at (decision) 
state i if we start at (non-decision) state j. It is known |13j that these probabilities 
can be calculated using the fundamental matrix, F = (I — Rq)^ 1 . More precisely, 



= QoF = Qo(I - RoY 



(3.17) 



and the computation requires a matrix inversion and a matrix multiplication. If we 
use classical methods, Q* can be calculated in 0(r 3 + r 2 s). Note that we could also 



apply more efficient algorithms, such as the Coppersmith- Winograd method. 
Q* the transition matrix of policy fi in the reduced problem should be 



Using 



Pa = Qa + Q*R, 



(3.18) 



This matrix encodes the idea that if we arrive at a non-decision state, we simulate 
the progress of the system until we arrive at a decision state. Fortunately, we do 
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not have to compute it for all possible policies, we only need to define the transition- 
probabilities accordingly. In order to achieve this, we can use 

T 

p(j\h u ) = P(j \ h u ) + ^v{k\i,u)Q* k (3.19) 
fc=i 

for all i,j>r and u £ U(i). Note that i and j are decision states. Computing the 
new transition-probability function needs 0(s 2 rm) operations. 

We should also modify the immediate-cost function, in order to include the ex- 
pected costs of those stages that we spend in the non-decision states, as well. It 
is known that the fundamental matrix contains information about the expected ab- 
sorbing times. More precisely, Fjk is the expected time spent in (non-decision) state 
j before arriving at a (decision) state (absorption), if the process started in (non- 
decision) state k [13] . Therefore, we define the new immediate-cost function as 

r r 

g{i,u) = g(i,u) + J2p(k\i,u)J2 F 3k9U), (3-20) 

k=l ] = 1 

for all i > r and u € U{i) 1 where we do not denoted the dependence of the cost- 
function on the actions for non-decision sates, since there is only one available action 
in each such state. Thus, g(j) = g(j,u), where u denotes the only available control 
action in state j. It is easy to see that computing the new cost- function needs 0{r 2 sm) 
operations, naturally, after we have calculated the fundamental matrix. 

We have only removed non-decision states, in which there is only one allowed 
action, consequently, it is trivial to extend a policy of the reduced problem to a policy 
of the original one, and there is a bijection between such policies. Since we defined the 
transition-probabilities and the immediate-cost function in a way that it mimics the 
behavior of the original problem, solving the reduced problem is equivalent to solving 
the original one. Summing all computations together, we can conclude that the time 
complexity of constructing the new, reduced SSP problem is 0(r 3 + r 2 sm + s 2 rm). 

3.9.2. Reducing the SSP formulation of Max-PageRank. Applying this 
result to the last, refined SSP formulation of the Max-PageRank problem (which also 
included damping and personalization), we can reduce the number of states of the 
problem to d + 1 by constructing another SSP problem as demonstrated above. The 
construction takes 0(n 3 + d 2 n + n 2 d) operations, after which the problem can be 
solved in 0(d 3 L) by using standard LP solvers [TT]. Consequently, we have 

Lemma 3.2. Max-PageRank can be solved in 0(n 3 +d 3 L) operations. 

4. PageRank Iteration. In this section we will provide an alternative solution 
to the Max-PageRank problem. We will define a simple iterative algorithm that in 
each iteration updates the configuration of the fragile links in a greedy way. Yet, as 
we will see, this method is efficient in many sense. For simplicity, we will only consider 
the case without damping (c = 0) and we will apply a new assumption as follows. 

(ABR) Bounded Reachability Assumption: We assume that the target node, v, can 
be reached from all nodes via a bounded length path of fixed edges. In other 
words, there is a universal constant k such that node v can reached from all 
nodes of the graph by taking at most k fixed edges. Constant k is universal 
which means that it does not depend on the particular problem instance. 
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The algorithm works as follows. It starts with a configuration in which each 
fragile link is activated. In iteration k it computes the expected first hitting time of 
node v if we start in node i and use the current configuration, that is 

H k (i) =E[min{i > 1 : X t = v } \X =i], (4.1) 

for all node i, where the transition matrix of the Markov chain (Xq, X\, . . . ) is P k 
defined by equation (|2.1[) using the adjacency matrix corresponding to the fragile link 
configuration in iteration k. Then, the configuration is updated in a simple, greedy 
way: a fragile link from node i to node j is activated if and only if H k (i) > H k (j) + 1. 
The algorithm terminates if the configuration cannot be improved by this way. We 
call this algorithm PageRank Iteration (PRI) and it can be summarized as follows: 



The PageRank Iteration Algorithm 



Input: A digraph Q = (V, £), a node v £ V and a set of fragile links T C £. 



1. 


k := 




% initialize the iteration counter 


2. 


F :=T 




% initialize the starting configuration 


3. 


Repeat 




% iterative evaluation and improvement 


4. 




= 1 T (/-Q,)- 1 


% compute the mean hitting times to v 


5. 


-Ffc+i 


:={{i,j)eT: 


H k (i) > H k {j) + 1} % improve the configuration 


G. 


k := 


k + l 


% increase the iteration counter 


7. 


Until F k - 


-x + F k 


% until no more improvements are possible 



Output: 1/H k (v), the Max-PageRank of v, and F k , an optimal configuration. 



Note that the expected first hitting times can be calculated by a system of linear 
equations [15 . In our case, the vector of hitting times, H k , can be computed as 

H k = t T (I-Q k )-\ (4.2) 

where Q k is obtained from P k by setting to zero the row corresponding to node v, 
namely, Q k = diag (1 — e v ) P k , where e v is the v-th n dimensional canonical basis 
vector. To see why this is true, recall the trick of Section [3~^l when we split node v into 
a starting node and an absorbing target node. Then, the expected hitting times of 
the target state can be calculated by the fundamental matrix [13] . If v can be reached 
from all nodes, then I — Q k is guaranteed to be invertible. Note that H k (v) = <p k (v), 
where tp k {v) is the expected first return time to v under the configuration in iteration 
k, therefore, the PageRank of v in the fc-th iteration is Tz k (v) = 1/H k (v). 

Theorem 4.1. PageRank Iteration has the following properties: 

(I) Assuming (AD) and (AR), it terminates finitely with an optimal solution. 
(II) Assuming (ABR), it finds an optimal solution in polynomial time. 

Proof. Part I: We can notice that this algorithm is almost the policy iteration (PI) 
method, in case we apply a formulation similar to the previously presented simple SSP 
formulation. However, it does not check every possible action in each state. It opti- 
mizes each fragile link separately, but as the refined SSP formulation demonstrates, 
we are allowed to do so. Consequently, PRI is the policy iteration algorithm of the 
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refined SSP formulation. However, by exploiting the special structure of the auxiliary 
states corresponding to the fragile links, we do not have to include them explicitly. 
For all allowed policy y, (for all configurations of fragile links) we have 



for all auxiliary states corresponding to a fragile link. Therefore, we do not have 
to store the cost-to-go of these states, because they can be calculated if needed. 

Notice that J^ k (i) = Hk(i), where /j,k is the policy corresponding to the con- 
figuration in iteration k. Thus, calculating is the policy evaluation step of PI, 
while computing Pfe+i is the policy improvement step. Since PRI is a PI algorithm, it 
follows that it always terminates finitely and finds an optimal solution [5] if we start 
with a proper policy and under assumptions (Al) and (A2). Recall that the initial 
policy is defined by the full configuration, Fq = T and that we assumed (AR), that 
is node v can be reached from all nodes for at least one configuration which means 
that the corresponding policy is proper. If this holds for an arbitrary configuration, it 
must also hold for the full configuration, therefore, the initial policy is always proper 
under (AR). Assumption (Al) immediately follows from (AR) and assumption (A2) 
follows from the fact that if the policy is improper, we must take infinitely often fixed 
or activated fragile links with probability one. Taking each one of these edges has 
unit cost, consequently, the total cost is infinite for at last one state. 

Part II: First, note that assumption (ABR) implies (AR) and (AD), therefore, 
we know from Part I that PRI terminates in finite steps with an optimal solution. 
Calculating the mean first hitting times, Hk, basically requires a matrix inversion, 
therefore, it can be done in 0(n 3 ). In order to update the configuration and obtain 
i*fc+i, we need to consider each fragile link individually, hence, it can be computed in 
0(d). Consequently, the problem of whether PRI runs in polynomial time depends 
only on the number of iterations required to reach an optimal configuration. 

Since we assumed (ABR), there is a universal constant k such that for all nodes 
of the graph there is a directed path of fixed edges from this node to node v which 
path has at most k edges. These paths contain fixed (not fragile) edges, therefore, 
even if all fragile links are deactivated, node v can still be reached with positive 
probability from all nodes. Consequently, all policies are proper (APP). It is easy to 
see that we can partition the state space to subsequent classes of states Si, . . . , S r , 
where r < k, by allocating node i to class S q if and only if the smallest length path 
of fixed edges that leads to node v has length q. This partition satisfies the required 
property described in Section 12.3.31 Because PRI is a PI variant, PRI terminates 
with an optimal solution in iterations bounded by a polynomial in L and rj~ 2K . Since 
77 = 1/m, where m < n is the maximum out-degree in the graph, r]~ 2K = 0(n 2K ), 
therefore, PRI runs in polynomial time under the Turing model of computation. □ 

5. PageRank Optimization with Constraints. In this section we are going 
to investigate a variant of the PageRank optimization problem in which there are 
mutual exclusive constraints between the fragile links. More precisely, we will consider 
the case in which we are given a set of fragile link pairs, C C T x T, that cannot be 
activated simultaneously. The resulting problem is summarized below. 

We saw that the Max-PageRank problem can be solved in polynomial time even if 
the damping constant and the personalization vector are part of the input. However, 
we will see that the Max-PageRank problem under exclusive constraints is already 




if Kfij) 



"off", 



on 



(4.3) 
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NP-hard, more precisely, we will show that the decision version of this problem is 
NP-complctc. In the decision version, one is given a real number p G (0, 1) and is 
asked whether there is a fragile link configuration such that the PagcRank is larger 
or equal to p. The proof is based on a reduction from the 3SAT problem. 



The Max-PageRank Problem under Exclusive Constraints 

Instance: A digraph Q = (V,£), a node i)6V,a set of controllable edges T C £ 

and a set C C T x T of those edge-pairs that cannot be activated together. 
A damping constant c £ (0, 1) and a stochastic personalization vector z. 

Task: Compute the maximum possible PageRank of v by activating edges in T 

and provide a configuration of edges in T for which the maximum is taken. 



Theorem 5.1. The decision version of the Max- PageRank Problem under 
Exclusive Constraints is NP-complete. 

Proof. The problem is in NP because given a solution (viz., a set of activated 
fragile links), it is easy to verify in polynomial time, e.g, by a matrix inversion, cf. 
equation (|2.5p . whether the corresponding PagcRank is larger or equal than p. 

We now reduce the 3SAT problem, whose NP-completeness is well known [10], to 
this problem. In an instance of the 3SAT problem, we are given a Boolean formula 
containing m disjunctive clauses of three literals that can be a variable or its negation, 
and one is asked whether there is a truth assignment to the variables so that the for- 
mula (or equivalently: each clause) is satisfied. Suppose now we are given an instance 
of 3SAT. We are going to construct an instance of the Max-PageRank Problem under 
exclusive constraints that solve this particular instance of 3SAT. 

We construct a graph having m + 2 nodes in the following way: we first put a 
node s and a node t. Figure it as a source node and a sink node respectively. Each 
clause in the given 3SAT instance can be written as j/j.i V yj.2 Vj/j.3, 1 < j < m, where 
Ujj is a variable or its negation. For each such clause, we add a node Vj between s 
and t, we put an edge from Vj to itself (a self- loop), we put an edge from s to vj, 
and we put three edges between Vj and t, labeled respectively with yj,i,yj,2, and yys- 
We finally add an edge from t to s. We now define the set of exclusive constrains, C, 
which concludes the reduction. For all pair {yj,i,yj> ,i>) such that yyi = yy ^ (i.e., 
is a variable and yy ,v is its negation, or conversely), we forbid the corresponding pair 
of edges. Also, for all pair of edges (yjj,yjj>) corresponding to a same clause node, 
we forbid the corresponding pair. This reduction is suitable, since our graph and the 
set C have a size which is polynomial in the size of the 3SAT instance. 

We claim that for c small enough, say c = l/(100m), it is possible to obtain an 
expected return time from t to itself which is smaller than 77 if and only if the instance 
of 3SAT is satisfiable. The reason for that is easy to understand with c = : if the 
instance is not satisfiable, there is a node Vj with no edge from it to t. In that case, 
the graph is not strongly connected, and the expected return time from t to itself 
is infinite. Now, if the instance is satisfiable, let us consider a particular satisfiable 
assignment. We first activate all edge which correspond to a literal which is true, 
and second, if necessary, we deactivate some edges so that for all clause node, there 
is exactly one leaving edge to t. This graph, which is clearly satisfiable, is strongly 
connected, and so the expected return time from t to itself is finite. 

Now if c 7^ is small enough, one can still show by continuity that the expected 
return time is much larger if some clause node does not have a leaving edge going to 
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node t. To see this, let us first suppose that the instance is not satisfiablc, and thus 
that a clause node (say, Vi), has no leaving edge. Then, for all I > 3, we describe a 
path of length I from t to itself: this path passes through s, and then remains during 
I — 2 steps in Vi, and then jumps to t (with a zapping). This path has probability 
(1 — c)tjt(1 — c)'~ 2 c. Therefore, the expected return time is larger than 



where we assumed that c = 1 / (100m) and the personalization vector is z = (1/n) 1. 
Note that we are allowed to determine c and z, since they are part of the input. 

Consider now a satisfiable instance, and build a corresponding graph so that for 
all clause node, there is exactly one leaving edge. It appears that the expected return 
time from t to itself satisfies E2 < 77.. To see this, one can aggregate all the clause 
nodes in one macro-node, and then define a Markov chain on three nodes that allows 
to derive a bound on the expected return time from vt to itself. This bound does not 
depend on m because one can approximate the probabilities m/(m + 2) and l/(m + 2) 
that occur in the auxiliary Markov chain by one so that the bound remains true. 
Then, by bounding c with 1/8 > l/(100m), one gets an upper bound on the expected 
return time. For the sake of conciseness, we skip the details of the calculations. 

To conclude, this is possible to find an edge assignment in the graph so that the 
pagerank is greater than p = 1/77 if and only if the instance is satisfiable. □ 

We have tried to keep the NP-hardness proof as short as possible. Several variants 
are possible. In the above construction, each clause node has three parallel edges 
linking it to the node t. This might seem not elegant, but it is not difficult to get rid 
of them by adding auxiliary nodes. Also, it is not difficult to get rid of self-loops by 
adding auxiliary nodes. Finally we have not tried to optimize the factor c = l/(100m), 
nor the bound on Ei- An interesting question is whether a reduction is possible if the 
damping factor c and the personalization vector z cannot depend on the instance. 

6. Conclusions. The task of ordering the nodes of a directed graph according 
to their "importance" arises in many applications. A promising and popular way to 
define such an ordering is to use the PageRank method [6] . The problem of optimizing 
the PageRank of a given node by changing some of the edges caused a lot of recent 
interest [TJ [12] . The paper considered a general problem of finding the extremal 
values of the PageRank a given node can have in the case we are allowed to control 
(activate or deactivate) some of the edges, which we referred to as "fragile links" . 

Our main contribution was that we proved that this general problem could be 
solved optimally in polynomial time under the Turing model of computation, even if 
the damping constant and the personalization vector were part of the input. The proof 
is based on reformulating the problem as a stochastic shortest path (SSP) problem (a 
special Markov decision process) and it results in a linear programming formulation 
that can be solved by standard techniques, e.g., interior point methods. 

Note that we do not need to assume that the graph is simple, namely, it can have 
multiple edges (and self-loops). Hence, our result can be generalized to PageRank 
optimization on weighted graphs, in case the weights are positive integers or rationals. 

Based on the observation that in some of the states of the reformulated SSP 
problem there is only one available action, we showed that the number of states (and 
therefore the needed computation to solve the problem) could be further reduced. 
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We also suggested an alternative solution, called the PagcRank Iteration (PPJ) 
algorithm, which had many appealing properties, for example, it always finds an 
optimal solution in finite steps and it runs in polynomial time, under the bounded 
reachability assumption. We left open the question of knowing whether PRI runs in 
polynomial time under milder conditions. The analysis and generalization of the PRI 
algorithm in case of damping and personalization is also a potential research direction. 

Finally, we showed that slight modifications of the problem, as for instance adding 
mutual exclusive constraints between the activation of several fragile links, might turn 
the problem NP-hard. We conjecture that several other slightly modified variants 
of the problem are also NP-hard, for example, the Max-PageRank problem with 
restrictions on the number of fragile links that can be simultaneously activated. 
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